Fokker-Planck and Langevin Equations from Forward— Backward Path Integral 



Hagen Kleinert* 
Institut fur Theoretische Physik, 
Freie Universitat Berlin, Arnimallee 14, 14195 Berlin, Germany 

Starting from a forward-backward path integral of a point particle in a bath of harmonic os- 
cillators, we derive the Fokker-Planck and Langevin equations with and without inertia. Special 
emphasis is placed upon the correct operator order in the time evolution operator. The crucial step 
is the evaluation of a Jacobian with a retarded time derivative by analytic regularization. 
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I. INTRODUCTION 



In 1963, Feynman and Vernon |l| set up a path integral for the study of a point particle in a thermal bath of 
harmonic oscillators. If the particle has a mass M, moves in a potential V(x), and is coupled linearly to a large 
number of oscillators Xi(t) of mass Mi and frequency tti, the probability to run from the spacetime point x a t a to 
x b t b is given by the forward-backward path 



\{xbtb\x a t a )\ 2 = / Vx + Vx- exp 



M 



(4 



- V(x+) + V(x-) -x+J2 c * x + + X -Y. ClX - 



(1.1) 



where Ci are coupling strengths. The bath oscillators are supposed to be in thermal equilibrium at a temperature T. 
This is taken into account by forming the thermal average of the bath oscillators. For a single oscillator, this is done 
by considering CiX± as external currents j± coupled to X±, and calculating the Gaussian integral 

(1.2) 



(1.3) 



Zo\j+,3-} = J dX b dX a (X b hf3\X a 0) n {X b t b \X a t a )£ (Xbtb\X a ta)o_ 
where (Xb h(3\X a 0)ci is the imaginary-time amplitude 

1 MQ 



(x b np\x a o) 



1 



y/2nh/M V sinh h/3 



exp 



2fi sinh h/3fl 



[{Xt + Xi) cosh h/3n - 2X b X a ] 



and (Xbtb\X a t a y n the path integral over the bath oscillator 



{X b tb\X a t a y Q = / VX{t)e*v\ - 



dt 



f { x*-tfxi) 



with a total classical action 



A 



1 



MCI 



2 sin fl(tb — t a 
1 



• [{Xl + X 2 a ) cos n(t 6 - t a ) - 2X b X a ] 



sin fl(tb - t a ) J ta 
and the fluctuation factor 

1 r 



dt[X a smCl(t b -t) + X b smn(t-t a )]j(t), 



Fnj{t b , t a ) 



y/2wih/M V sin fl(t b - t a ) 
x exp 



dt dt' S mCl(t b -t)sma(t'-t a )j(t)j(t') 
nMUsmil(t b — t a ) J t J t 



(1.4) 



(1.5) 



(1.6) 
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The result of the thermal average in Eq. (1.2) is 

Zo\j+,j-] = exp | - ^ J dt J dt'e(t - t') 

x "(j+ - 3-)(t)G R (t,t>)( J+ + j_)(t') + (j + -j-)(t)A(t,t')(j + -j-)(t') , 



(1.7) 



where A(t, t') and G R (t, t') are the expectation values {{X{t), X(t')}) and Q{t— t')([X(t), X(t')]) in operator language. 
The are real and imaginary parts of the time-ordered Green function and G(t,t') — <d(t — t'){TX(t), X{t')): 



n 



n cosh- [h(3 - i(t - t')] 
G(t, t') = A(t, t') + G R (t, t') = 2 



2MCl 



sinh 



t > t'.. 



(1.8) 



which is the analytic continuation of the p erio dic imaginary-time Green function to r = it. 

The thermal average of the probability ([bl]) is then given by the forward-backward path integral 



\(Xbt b \Xata)\ 2 = / Vx+{t) / VX-{t) 



x exp < — 



dt 



M 

-{xl-x 2 _)-{V{x + )-V{x.)) 



where expz„4 FV [a;+, x-]/h is the Feynman- Vernon influence functional defined by 
Zq[x+, x_] = exp {iA FV [x+, X-]/h} = exp {iAj^ [x+, x-]/h + iA^lx^, x-]/hj 



(1.9) 



(1.10) 



exp 



2Ti 



dt / dt'e(t-t') {x + -x-)(t)G£{t,t')(x+ +x-)(t') + {x+ - x-){t)A h (t,t'){x+ -x_)(i') 



with G^(t,t') and A\,{t, t') being commutator and anticommutator functions of the bath. The first and second parts 
of the exponents have been distinguished as dissipative and fluctuating parts A F J / [x + , x_] and „4 FV [x + , x_] of of the 
effective influence action A FY [x + , xJ\. 

The bath functions G£(t,t') and A h (t,t') are sums of correlation functions of the individual oscillators of mass Mi 
frequency Oj, each contributing with a weight cf . Thus we may write 

/OO 7 / 

^* h (u,')i S mu'(t-t'), 
- oo 27T 

A h (t,t') =^<$({Mt),X<V)})p = fi J |%V)coth 
where <Tb(u>') is the spectral density of the bath 



hoj' 

2k B T 



COS0j'(t — t'), 



(1.11) 



* h (u,') EE 27T gjjfe^' - ^ - ^ + ^ 



For a discussion of the properties of the influence functional, we introduce an auxiliary retarded function 

1 r°° dw£(w) -iw(t-f) 



j(t-t') = e(t-t') 



M 



2ir u> 



and write 



with 



MAuj 2 



G R (t, t') = ihMj(t - t') + i Auj 2 S(t - t 1 ) 

dui' CTb(w') 



2tt w' 



E 



(1.12) 



(1.13) 



(1.14) 



(1.15) 
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Inserting the first term in the decomposition (1.14) into (1.11), the dissipative part of the influence functional can be 
integrated by parts in if and becomes 



Ajy[x+,X-} = — - / ill I dl' [.>■--. v )(/ )-:(!'- I )!./•+ + r II/') 
M 

+ — / dt{x+-x-)(t)j(t b -t)(x + +x-)(t a ). 



The 5-function in (1.14) contributes to .Ap V [x + , £_] a term 



AA\ oc [x+, x-] 



M 



dtAuj 2 {x 2 + - xi)(t), 



which may simply be absorbed into the potential terms in the path integral (l.S) renormalizing them to 



dt [V ron {x + ) - V len (x-)] 



(1.16) 



(1.17) 



(1.18) 



The odd bath function (Tb(w') can be expanded in a power series with only odd powers of u>'. The lowest approxi- 
mation 



CTb(w') « 2M7W', 



(1.19) 



describes Ohmic dissipation with some friction constant 7. For frequencies much larger than the atomic relaxation 
rates, the friction goes to zero. This behavior is modeled by the Drude form of the spectral function 



a{w') « 2M7W 



I _fD_ 

0J 2 n + bJ r 



Inserting this into Eq. ( 1.13Q , we obtain the Drude form of the function "f(t)\ 

7 g(i) = e(t) 7 ^e- D l*l. 
The superscript emphasizes the retarded nature. This can also be written as a Fourier integral 

J-ioo 27T 

with the Fourier components 

iui D 



7^0=7-7 



u>' + ilud 



(1.20) 



(1.21) 



(1.22) 



(1.23) 



The position of the p ole in the lower half-plane ensures the retarded nature of the friction term by producing the 
Heaviside function in (1.21). 

In the Ohmic limit (1.19), the dissipative part of the influence functional simplifies. Then 7(4) becomes narrowly 
peaked at positive t may be expressed in terms of a left-sided ^-function as 

7 (i)-27**(t), (1.24) 

whose superscript R indicates the retarded asymmetry of the ^-function, which has the property that 

dt@(t)S R {t) = 1. (1.25) 



With this, (1.16) becomes a local action 



Ajj , x_ 



M 



7 / dt{x+ - x_)(x + + i_) - -7T7(x+ - xt)(t a ) 



(1.26) 
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The right-sided nature of the 5-function causes an infinitesimal negative shift in the time argument of the velocities 
(x + +X-)(t) with respect to the factor (x + —X-)(t), indicated by the superscript R. It expresses the causality of the 
friction forces and will be seen to be crucial in producing a probability conserving time evolution of the probability 
distribution. 

The second term changes only the curvature of the effective potential at the initial time, and can be ignored. In 
t he fi rst term it is important to observe that the retarded nature of the dissipative term and of the function j(t) in 
(1.13) ensures that the velocity term (i+ +x-)(t) lies before (x+ — X-)(t) in a time-sliced path integral. This ensures 
the causality of the friction forces. 

It is useful to incorporate the slope information (1.19) also into the bath correlation function Ab(i, t') in (1.11), 
and factorize it as 



where 



A h (t,t') = 2Mjk B TK(t,t'), 
K(t,t') = K(t-t') = 



(1.27) 



2Mjk B T 

r duj ' k( <\ 

loo ^ K[UJ ^ 



iu}' (t— t') 



with Fourier transform 



1 u(lo') Huj' , tuv' 
■ coth ■ 



2M 7 u' 2k B T 
which in the limit of a purely Ohmic dissipation simplifies to 



K(fi/) = K° hm (uj') - 



coth 



2k B T' 



2k B T 2k B T 

The function K(u)') has the normalization K(Q) = 1, giving K(t — t') a unit temporal area: 



dtK(t-t') = 1. 



In the classical limit fi — » 0, 



and 



K{u') = 



K(t - t') = 



2oj, 



,-uin{t-t') 



(1.28) 



(1.29) 



(1.30) 



(1.31) 



(1.32) 



(1.33) 



In the limit of Ohmic dissipation, this becomes a <5-function. Thus K(t — t') may be viewed as a J-function broadened 

by quantum fluctuations and relaxation effects. 

With the function K(t,t'), the fluctuation part of the influence functional in fll.10 ), (1.11), (1.9) becomes 



.M~fk B T 



dt 



dt' (x+ - x-){t)K(t,t') (x+ - x-)(t'). 



(1.34) 



Here we have used the symmetry of the function K(t, t') to remove the Heaviside function Q(t — t') from the integrand, 
extending the range of i'-integration to the entire interval (t a , tb). 

In the Ohmic limit, the probability of the particle to move from x a t a to Xbtb is given by the path integral 



|(a:&f&|a; t )| 



Vx + (t) I Vx-{t) 
M 



i 

x exp . 

n 

x exp < —i 



dt 



(x\-xt)-{V{x + )-V{x_)) 



M7 
~2 

M"fk B T 



dt^-(x + -X-)(t)(x + + i-) R (t) 



dt 



dt 1 (x + -x-)(t) K(t,t') (x + -x_){t') 



(1.35) 



4 



This is the closed-time path integral of a particle in contact with a thermal reservoir. 

The paths x+(t ), X—it) may also be associated with a forward and a backward movement of the particle in time. For 
this reason, ( 1.35[ ) is also called & forward-backward path integral. The hyphen is pronounced as minus, to emphasize 
the opposite signs in the partial actions. 

It is now convenient to change integration variables and go over to average and relative coordinates of the two paths 
x+, X-\ 



X 



(x++x.)/2, 
x + — X-. 



Then 



1.35|) becomes 



\{x b t b \x a t a )Y 



Vx{t) I Vy(t) exp 



dt 



M {-yx + 1V x R ) + V (x + |) - V (x - |) 



Mjk B T 



dt 



dt'y(t)K(t,t')y(t') 



(1.36) 



(1.37) 



II. FOKKER-PLANCK EQUATION 



At high-temperatures, the Fourier transform of the K ernel K(t, t') in Eq. (1.30) tends to unity such that K(t,t') 
becomes a ^-function, and the bath correlation function (1.27) becomes approximately 



Ab(M') « ~w6(t-t'), 
where we have introduced the constant proportional to the temperature: 

w = 2Mjk B T, 
which is related to the so-called diffusion constant 

D = k B T/M-y 

by 

w = 2 1 2 M 2 D/T. 



(2.1) 

(2.2) 
(2.3) 
(2.4) 



Then the path integral (1.37) for the probability distribution of a particle coupled to a thermal bath simplifies to 



P(x b t b \x a t a ) = \(x b t b \x a t a )\ 



Vx{t) / Vy{t) 



x exp 



dt y[Mx + M-f± R + V'(x)] 



w 

2T? 



dty< 



(2.5) 



The superscript R records the infinitesimal backward shift of the time argument as in Eq. ( 1.26| ). The y- variable can 
be integrated out, and we obtain 



P{x b t b \x a t a ) 



Vx(t) exp 



1 

2w 



dt [Mi + Mjx R + V'(x)f 



This looks like a cuclidcan path integral associated with the Lagrangian 

1 



2w 



[Mx + M-fi + V'(x)} 2 



(2.6) 

(2.7) 
. The 



The solution of such path integrals with squares of second time derivatives in the Lagrangian is given in Ref. 
result will, however, be different, due to time-ordering of the x^-term. 

Apart from this, the Lagrangian is not of the conventional type since it involves a second time derivative. The 
action principle 8 A — now yields the Euler-Lagrange equation 

d dL d 2 dL „ 



dL 

dx 



dt dx dt 2 d'x 



0. 



This equation can also be derived via the usual Lagrange formalism by considering x and x as independent generalized 
coordinates x, v. 
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III. CANONICAL PATH INTEGRAL FOR PROBABILITY DISTRIBUTION 



It is well-known that a path integral satisfies a Schrodinger type of equation. For the path integral (2.6) this is 
known as a Fokker- Planck equation. The relation is established (see the textbook by rewriting the path integral 
in canonical form. Treating v = x as an independent dynamical variable, the canonical momenta of x and v are (see 
Section 17.3 of the textbook Ref. {§) 

V = i^r = i—[Mx + M^yx + V'(x)} 
ox w 

= i — -[Mv + Mvv + V'(x)], 
w 

Pv = i-^r = -p- (3.1) 

ox 7 



The Hamiltonian is given by the Legendre transform 



H(p,p v ,x,v) = L e (x,x) — y^g^^t (3-2) 

i—i % 

= L e (v, v) + ipv + ip v v, (3-3) 



where v has to be eliminated in favor of p v using (3.1). This leads to 



H(p,p V) x,v) = -^TjiPl -iPv[jv+ jj v '( x )} +ipv- (3.4) 



The the canonical path integral representation for the probability reads therefore 
P(x b t b \x a t a ) = [vx [vv 1 Dl ' r 



2tt J J 2vr 

xexpj^ dt [i(px +p v v) - H(p,p v ,x,v)] j. (3.5) 

It is easy to verify that t he p ath integral over p enforces v = x = x, after which the path integral over p v leads back 
to the initial expression ( |2.6| ). We may keep the auxiliary variable v(t) as an independent fluctuating quantity in all 
formulas and decompose the probability P(xbtb\x a t a ) with respect to the content of v as an integral 

/oo />oo 
dvb / dv a P(xbVbtb\x a V a t a )- (3.6) 
-oo J —oo 

The more detailed probability on the right-hand side has the path integral representation 

Vp v 



P{x b Vbt b \x a v a t a ) = \(xbVbtb\x a v a t a )\ 2 = J Vx J " J Vv J 



2tt 



xexp|^ dt [i(px+p v v) — H(p,p v ,x, , (3.7) 

where the end points of v are now kept fixed at Vb — u(tfe), v a = v(t a ). 

We now use the relat ion between a canonical path integral and the Schrodinger equation to conclude that the 
probability distribution (3.7) satisfies the Schrodinger-like differential equation:]] 



H(p,p v ,x,v)P(xvtb\x a v a ta) = -dtP(xvt\x a v a ta)- (3.8) 
This is the Fokker-Planck equation in the presence of inertial forces. 



1 See the review paper by S. Chandrasekhar, Rev. Mod. Phys. 1,5, 1 (1943). 
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At this place we no te t hat when going over from the classical Hamiltonian (3.4) to the Ha milto nian operator in the 
differential equation (|3.8| ) there is an operator ordering problem. When writing down Eq. (3.8) we do not know in 
wh ich order the momentum p v must stand with respect to v. If we were dealing with an ordinary functional integral 
in ( |2.6|) we would know the order. It would be found as in the case of the electromagnetic interaction to by symmetric: 
-(p v v + vp v )/2. 

On physical grounds, it is easy to guess the correct order The differential equation ( |3.8| ) has to conserve the total 
probability J dx dvP(x v tb\x a v a t a ) for all times t. This is guaranteed if all mome nta s tand to the left of all coordinates 
in the Hamiltonian operator. Indeed, integrating the Fokker-Planck equation (3.8) over x and v, only a left-hand 
position of the momentum operators leads to a vanishing integral, and thus to a time independent total p roba bility. 
We suspect that this order must be derivable from the retarded nature of the velocity in the term yx R in (2.5). The 
proof that this is so is the essential point of this paper, by which it goes beyond an earlier treatment of this subject 
in Ref. §. 



IV. SOLVING THE OPERATOR ORDERING PROBLEM 



Since the ordering problem in the Hamiltonian operator associated with (3.4) does not involve the potential V(x), 
we study this problem most simply by considering the free Hamiltonian 



H (p,Pv,x,v) 



2M 2 



P v — ijPvV + ipv. 



(4.1) 



which is associated with the Lagrangian path integral 



Pa{xbtb\x a ta) oc I Vx(t) exp \ -— J dt [Mx + M-f± R } 2 



(4.2) 



We furthermore may concentrate on the probability with x b = x a = 0, and assume t b — t a to be very large. Then the 
frequencies of all Fourier decompositions are continuous. 

Forgetting for a moment the retarded nature of the velocity x, the Gaussian path integral can immediately be done 
and yields 



P o {0t b \0t a ) ocDet- 1 (-9 t 2 - 7 9 t ) 



oc exp 



-{tb-ta) I P l0g(w' 2 - i 7 w') 



(4.3) 



The integral can be evaluated in analytic regularization according to the rule 

^log(^ 7 ) = ^, 



(4.4) 



which follows directly by integrating in 7 the sym metric decomposition of the integral J dio' /(uj' 2 + 7 2 ) = 7r/| 7 |- Hence 
we obtain for the functional determinant in ( |f ,3| ): 



Det (-dt - 7<9 t ) = Det (id^Bet (id t - 17) = exp [Tr log(«<9 t ) + Tr log(id t - i 7 )] 



exp 



(h-Q- 



(4.5) 



and thus 



P {0t b \0t a ) ocexp -(t b -t a ) 



(4.6) 



This corresponds to an energy 7/2 and an ordering —i"/(p v v + vp v )/2 in the Hamilton operator. 

We now take the retarded time argument of x R into account. Spe cifically, we replace the term ^yx R in (O) by 



[ dtd t' j/7n(f — t') x(t) containing the retarded Drude function (1.21) of the friction. The the frequency integral in 
([4.3|) becomes 



dio 
2^ 



log ( uj 12 - 7-= 



div 
2^ 



-\og(uj' + iuj D ) + log (lo 12 + iuj d -juj'u!d)] , 



(4.7) 
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where we have omitted a vanishing integral over logo/ by (4.4). We now use (4.4) to find 

\ud\ 



dio 
2^ 



- log(a/ + iuod) + log {yj' 2 + iu'ujd - 7^z?)] = 



2 



1^2 



0. 



where — iu>i2 are the solutions of the quadratic equation 

uj' 2 + ioj£) — juj'llJe) 
For a large Drude frequency lud, they are given by 

ojx = cod — 7, ^2 = 7 



(4.8) 



(4.9) 



(4.10) 



Inserting these into ( |4.8| ) we find a vanishing integral rather than 7 in (4J3), and thus a functional determinant 

Det {-d 2 t - 7 0« ) = exp [Tr log(-d 2 - 7 df )] = 1, (4.11) 

instead of (|4.5|). The notation symbolizes now the specific retarded functional matrix a la Drude with a large 
lu d : 



With the determinant (4.11), the probability becomes a constant 

Po(0< 6 |0t a ) = const, 



(4.12) 



(4.13) 



This shows that the retarded nature of the friction force has subtracted an energy 7/2 from the energy in ( |l.(j ). 
With the ordinary path integral corresponding to a Hamilton operator with a symmetrized term — i(p v v + vp v )/2, the 
subtraction of 7/2 has changed this to —i"fp v v. 



Note that the opposite case of an adva nced velocity term x in (4.2) would be approx ima ted by a Drude function 
7o(t) which looks just like 7o(t) in (1.23), but with negative ujo- The right-hand side of (4.8) becomes now 27 rather 
than zero, The corresponding formula for the functional determinant is 



Det {-d 2 - jdf) = exp [Tr log(-d t 2 - 7^)] = exp [(t b - t a ) 7 ] 



(4.14) 



where ^ydf stands for the advanced version of the functional matrix (4.12) in which u>d is replaced by —u>d- Thus we 
find 

P (flt b \Ot a )ocexp[-(t b -ta)i\, (4.15) 

with an additional energy 7/2 with respect to the ordinary formula ( [4.6[ ). This corresponds to the opposite (unphysical) 
operator order —ijvp v in Ho, which would violate the probability conservation of time evolution twice as much as the 
symmetric order. 

The above formulas for the functional determinants can easily be extended to the slightly mor e general case where 
V(x) is the potential of a harmonic oscillator V(x) = Muj 2 x 2 /2. Then the path integral (2.6) for the probability 
becomes 



Po(x b t b \x a t a ) cx / Vx(t) exp 



1 

2w 



dt [Mx + M^x R + Uqx) 



(4.16) 



which we evaluate at x b — x a — 0, where it is given by the properly retarded expression 

P (0 1 b 1 t a ) oc Det - 1 ( - d 2 - jd t + lo 2 ) 

du 



oc exp 



■(fib-ta) 



2tt 



log(w - i-fu>' - w 2 ) 



The roots of the argument of the logarithm lie now at —10^12 with 



W12 




(4.17) 



(4.18) 
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Using the a naly tically regularized formula (4.4), the integral in (4.17) yields simply 7, and we obtain the generalization 
of formula ( 15 ) 



Det {-d 2 - jd t - u 2 ) = exp [Tr log(-d t 2 - jd t - wg) 



exp 



{h-ta)~ 



(4.19) 



The frequency luq has no influence upon the result. This formula can be generalized further to time-dependent 
coefficients 



Det [~d 2 ~ j(t)dt - Q 2 (t)] = exp {Tr log [-d 2 - -yd t - ft 2 (i)] } = exp 

This follows from the factorization 

Det [-d 2 - j{t)d t - fl 2 (t)] = Det [d t + Det [d t + fl 2 (t)} 

with 

Qi(t) + n 2 (t) = 7 (0, d t n 2 {t) + ni(t)n a (t) = n 2 (t), 

using the formula 

"1 



Det [d t + j{t)} = exp 



2 Jt 



The probability of the general path integral (2.6) without retardation of the velocity term is therefore 



P o (0t b \0t a ) cxexp 



~(tb-tah 



(4.20) 

(4.21) 
(4.22) 

(4.23) 
(4.24) 



as in (|4.6|). 

Let us now introduce retardation of the velocity term by using the (independent Pru de ex pression ( 1.23j ) for the 
friction coefficient. First we consider again the harmonic path integral (4.16), for which (4.17) becomes 



/^(0/,|()/„) x exp{ (//,-/„)/ ^ log [lu' 2 -i^{lo)lu' -u; 2 ] 



For a large Drude frequency u>n 3> 7 the roots are now 

7 



W12 



1± ll- 



4m„ 



7 



w 3 = uj d - 7. 



Using once more formula (4.4), we see that 7 and the mo-terms disappear, and we remain with 

P o (0t b \0t a ) = const. 

This implies a unit functional determinant 

Det (d t 2 -M 7 <9f + oj 2 q ) = 1, 



(4.25) 



(4.26) 



(4.27) 



(4.28) 



in contrast to the unretarded determinant (4.19). By analogy with (4.21), the also the general retarded determinant 
is independent of j(t) and fl(t). 



Det [_02_ 7 (t)0«_n 2 (t)] =1. 
In the advanced case, we would find similarly 



Det [-d 2 - j{t)df - n 2 (t)] = exp 



dt-f(t) 



(4.29) 



(4.30) 



All three determinants are correct also for finite time intervals, due to the solvability of the first-order differential 
equation by means of an integrating factor. 



Note that the retardation prescription can be avoided by a trivial additive change of the Lagrangian (2.7) to 



L e {x, ±) = -L [a + M 7 i + V'{xf - J. 
From this the path integral can be calculated with the usual time slicing. 



(4.31) 
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V. STRONG DAMPING 



For 7 3> V"(x)/M, the dynamics is dominated by dissipation, and the Lagrangian (2/7) takes a more conventional 
form in which only x and x appear: 



L e(a; ,i) = l-[M 7 i ; « + y'(x)] 2 = 2 L 



where x R lies slightly before V'(x(t)). The probability 



P(cc 6 t h |a; a i a ) 



Vx exp 



dtL e (x, x R ) 



(5.1) 



(5.2) 



looks like an ordinary euclidean path integral for the density matrix of a particle of mass M = 1/2D. As such it 
obeys a differential equation of the Schrodinger type. Forgetting a moment about the subtleties of the retardation, 
we introduce an auxiliary momentum integration and go over to the canonical representation of (5.2): 



P(x b t b \x a t a ) = I Vx I exp 



dt 



v 2 1 
ipx-2D^-+ip—V'(x) 
2 M7 



This probability distribution satisfies therefore the Schrodinger type of equation 

H(p b , x b )P(x b t b \x a t a ) = -dt b P{x b t b \x a t a ) 

with the Hamiltonian operator 

H(p,x)=2D^-ip^-V'(x). 
2 Mj 



(5.3) 



(5.4) 



(5.5) 



In order to conserve probability, the momentum operator has to stand to the left of the potential term. Only then 
does the integral over x b of Eq. (|5~4|) vanish. Equation (5.4) is the overdamped or ordinary . 

Without the retardation on x in ( |5.2| ), the path integral would certainly give a symmetrized operator —i[pV'(x) + 
V'{x)p]/2 in H. This follows from the fact that the coupling (l/2DM r y)xV (x) looks precisely like the coupling of 
a particle to a magnetic field with a "vector potential" A(x) = (l/2DM^f)V'(x). In this case we can also perform 
immediately the path integral (|5.2|) 



(5.6) 



P Q {x b t b \x a t a ) oc / Vx{t) exp <j - j- / dt [M 7 i + V'{x)] 2 



at x b = x a = 0, which is given by 



P (0t b \0t a ) oc Det - 1 [d t + V'(x)/M 1 ] , 



where from formula (4.23) 



Det [d t + V"(x)/M~f] = exp 



dtV"{x)/2Mj 



(5.7) 



(5.8) 



The effect of retardation of the velocity in (5.1) is obvious since the trivial retarded determinant (4.29) is independent 
of the strength of the damping: 



Det [flf + V"(x)/Mj] = 1. 

In the advanced case one would have 

Det [df + V"(x)/Mj] =exp J dtV"(x)/Mj 

For the differential equation (|]4|), the difference between the ordinary and the retarded results ( |5.8D and ( |5.9| ) 
implies that the initially symmetric operator order — i[pV'(x) + V'(x)p]/2 in H is changes into —i[pV'(x) + V'(x)p]/2 — 
V"(x)/2 — ipV'(x), as necessary for conservation of probability. 



(5.9) 



(5.10) 
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As in Eq. (4.31) we can avoid the retardation of the velocity by adding to the Lagrangian ( |5.l| ) a term containing 
the second derivative of the potential: 



L c (x, x) 



1 

AD 



1 



2Af 7 



V"(x). 



(5.11) 



From this the path integral can be calculated with the same slicing as for the gauge-invariant coupling of a magnetic 
vector potential (see Sections 10.6 and 11.3 in the textbook H) 



P (x b t b \x a t a ) oc / Vx(t) exp 



1 



AD 



1 



AD 



V'{x) 



-i 2 



M7 



V"{x) 
2M7 



(5.12) 



VI. LANGEVIN EQUATIONS 



For for high jT. the forward-backward path integral (1.37) has only small fluctuations of y, and K(t,t') becomes 
a (5-function. Then we can expand 



v ( x + 1) v ( x ~ I) ~ yV ' {x) + Ya v "' {x) + ■ ■ 

keeping only the first term. We further introduce an auxiliary quantity rj(t) by 

r){t) = Mx R {t) + M~/x(t) + V'{x(t)). 



(6.1) 



(6.2) 



With this, the exponential function in (1.37) becomes after a partial integration of the first term using the endpoint 
properties y(t b ) = y(t a ) = 0: 



exp 



w 

at yr\ ^ 

y 1 2h 2 



dty 2 {t)\, 



(6.3) 



where w is the constant (|2.2|). The variable y can obviously be integrated out and we find a probability distribution 



P[tj] ocexpj-— jf dtn\t)^, 
The expectation value of an arbitrary functional of F[x] can be calculated from the path integral 

(F[x])r,=Af [ Vx P[rj\F[x], 



(6.4) 



(6.5) 



where the normalization facto r M is fixed by the condition ( 1 ) = 1. By a change of integration variables from x(t) 
to r](t), the expectation value (6.5) can be rewritten as a functional integral 

(6.6) 



(6.7) 



(F[x]) n =Af J V V P[ V ] F[x], 

In principle, the integrand contains a factor J _1 [cc], where J[x] is the functional Jacobian 

J[x] = Det [Srj(t)/Sx(t')} = det[M5 t 2 + Mjd? + V"(x(t))]. 



However, in Eq. (1.2S ) we have seen that the retarded determinant is unity, thus justifying its omission in ( |6.6| ). 

The path integral (6.6) may be interpreted as an expectation value with respect to the solutions of a stochastic 
differential equation (6.2) driven by a Gaussian random noise variable rj(t) with a correlation function 



{v{tW)) T =wK{t-t'), 



(6.8) 



The stochastic differential equation ( |6.2[ ) has thus been derived from the forward-backward path integral at high 
temperatures. 
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In the overdamped limit, the classical Langevin equation with inertia (6.2) reduced to the overdamped Langevin 
equation-n 



T)(t) = M~f±{t) + V'{x{t)). 



(6.9) 



The noise average is formed as in (p. 



VII. SUPERSYMMETRY 



Recalling the origin ( |5.8p of the extra last term in the exponent of the path integral (5.12), this can be rewritten in 
a slightly more implicit but useful way as 



P (x b t b \x a t a ) cc / Vx(t)Bet 



V"(x) 



1 



exp < — — — / dt —- 



AD 



1 



AD 



V'(x) 



(7.1) 



In the form (7.1), the time ordering of the velocity term is arbitrary. It may be quantum mechanical, but equally well 
retarded or advanced, as long as it appears in the same way in both the Lagrangian and determinant. An interesting 
structural observation is possible by generating the determinant with the help of an auxiliary fermion field c(t) from 
a path integral over c(t): 



dBb[dt + V"(x{t))/Mi\ oc / VcVce- I ' dtm[M-yd t +v" (x(t))]c(t) 



(7.2) 



In quantum fi eld theory, such auxiliary fermionic fields are referred to as ghost fields. With these we can rewrite the 
path integral (5.2) for the probability as an ordinary path integral 



P(x b t b \x a t a ) — J Vx J VcVc exp{-Aps[x,c,c}} . 



where ^4ps is the euclidean action 

1 



A 



PS 



2L>MV 



dt { \ [Myx + V'(x)f + c(t) [Mjd t + V"{x(t))} c(t) 



(7.3) 



(7.4) 



first written down by Parisi and Sourlas.f] This action has a particular property: If we denote the expression in the 
first brackets by 



U x = Mjd t x + V'(x), 



the operator between the Grassmann variables in (7.4) is simply the functional derivative of U c 



su 

U xy = -^-=Mjd t + V"(x). 



Thus we may write 



ApS= 2Dj t dt 



-U 2 x +c(t)U xy c(t) 



(7.5) 



(7.6) 



(7.7) 



where U xy c(t) is the usual short notation for the functional matrix multiplication J dt'U xy (t,t')c(t'). The relation 
between the two terms makes this action super symmetric. It is invariant under transformations which mix the Fermi 



2 The stochastic movement is now a so-called Wiener process. 

3 G. Parisi and N. Sourlas, Phys. Rev. Lett. 43, 744 (1979); Nucl. Phys. B 206, 321 (1982). 
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and Bose degrees of freedom. Denoting by e a small anticommuting Grassmann variable, the action is invariant under 
the field transformations 

Sx(t) = ec(t) + c(t)e, (7.8) 

5c(t) = -eU x , (7.9) 

Sc(t) = U x e. (7.10) 

The invariance follows immediately after observing that 

SU X = eU xy c{t) + c(t)U xv e. (7.11) 

Formally, a similar construction is also possibl e for a particle with inertia in the path integral ( p.6| ), which is an 
ordinary path integral involving the Lagrangian (4.31). Here we can write 



1 



P(x b t b \x a t a ) =M J Vx(t)J[x]exp^-^- j dt[Mx + M~/x + V'(x)] 2 j. (7.12) 

where J[x] abbreviates the determinant 

J[x] = det[Md 2 + Mjd t + V"(x(t))} (7.13) 

which is known from formula ( 4.2CQ . The path integral ( 7.12j ) is valid for any ordering of the velocity term, as long 
as it is the same in the exponent and the functional determinant. 

We may now express the functional determinant as a path integral over fermionic ghost fields 



J[x] = det[Md? + M 7 d t + V"{x(t))] oc / PcDce"/ dt5(f) [ Ma ' +M79i+1/ " (l(f)) l cW I (7.14) 
and rewrite the probability P(x b t b \x a t a ) as an ordinary path integral 

P(x b t b \x a t a ) (x J Vx(t) J VcVcexp{-A KS [x 1% c}}, (7.15) 



where A[x, ,c] is the euclidean action 



^ ks [.t 71 c] = f " dt(^[Mx+Mjx+V'(x)f +c(t) [Mdf + Mjdt + V'ixit))] c(t) 



(7.16) 



This formal expression contains subtleties arising from the boundary conditions when calculating the Jacobian (7.14) 
from the functional integral on the right-hand side. It is necessary to fa ctori ze the second-order operator in the 
functional determinant and express each factor as a determinant as in (7.14). At the end, the action is again 
supersymmetric, but there are twice as many auxiliary Fermi fields ||. 



VIII. STOCHASTIC QUANTUM LIOUVILLE EQUATION 



At lower temperatures, where quantum fluctuations become important, the forward-backward path integral (1.37) 
does not allow us to derive a Schrodinger-like differential equati on fo r the probability distribution P(x v t\x a v a tt). To 
see the obstacle, we go over to the canonical representation of (1.37): 

\(x b t b \x a t a )\ 2 = J Vx{t) Vy{t) J exp jf " dt [px + p y y - H T ] 

(8.1) 

where 

1 w 
H T = jjPyPx + IPyV + V{x + y/2) - V{x - y/2) - i^yKy (8.2) 
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plays the role of a Hamiltonian. Here Ky(t) denotes the product of the functional matrix K(t,t') with the functional 
vector y(t') defined by Ky(t) = J dt' K(t,t')y(t'). After omitting the y-integrations at the endpoints, we obtain a 
path integral representation for the product of amplitudes 



P(xbybtb\x a y a t a ) = (xb + yb/2 h\x a + y a /2 t a )(x b - y b /2 t b \x a - y a /2 t a )* . 
The Fourier transform of the diagonal elements of this with respect to y is the Wigner function 



1 f°° 

W(x,p,t) = — J e?™»' n p{xyt\xyO). 



•3) 



(8.4) 



When considering the change of p(x y t\x a y a t a ) over a small time interval e, the momentum variables p and p y have 
the same effect as differential operators —id Xb and —id Vb , respectively. The last term in Ht, however, is nonlocal in 
time, thus preventing a derivation of a Schrodinger-like differential equation. 

The locality problem can be removed by introducing a noise variable rj(t) with the correlations function 



(v(tW))r, = ^K- 1 (t,t'). 

Then we can define a temporally local ^-dependent Hamilton operator 

H n = -^(p x + 7 y) p y + V(x + y/2) - V{x - y/2) - y V 



.5) 



(8.6) 



which governs the evolution of ry-dependent versions of the amplitude products (8J5) via the stochastic Schrddinger 
equation 



ih~d t p v (xyt\x a y a t a ) = H v p v (xyt\x a y a t a ). 



•7) 



Averaging this equation over r\ using ( 3J) ) yields for y a = y b = the same probability distribution as the forward- 
backward path integral (1.37): 

\{x b t b \x a t a )\ 2 = p(x b 0t b \x a 0t a ) = (p(x b 0t b \x a y a t a )} v (8.8) 



At high temperatures, the noise averaged stochastic Schrodinger equation (B.7) takes the form 

iWtp{xyt\x a y a ta) = Hp(xyt\x a y a t a ), 



(8.9) 



where H is the Hamiltonian associated with the Lagrangian in the forward-backward path integral (1.37): 

H = jjPyPx + lypy + V(x + y/2) - V{x - y/2) - i^y 2 . 



(8.10) 



In terms of the separate path positions x± = x ± y/2 where p x — 9+ + 9_ and p y = (9+ — 9_)/2, this takes the more 
familiar form 



1 



H = ^ {Pi - pi) + V(x+) - V(x-) + \{x+- x-){p+ - P-) - ihk(x+ - x_f. 
In the last term we have introduced a useful quantity, the so-called decoherence rate per square distance 



A 



2Ti Ti 

It is composed of the damping rate 7 and the squared thermal length 



as 



A : 



2irTi 2 (3/M 

27T7 



.11) 



(8.12) 



.13) 



(8.14) 
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and controls the rate of decay of interference peaks.^ 

Note that the order of the operators in the mixed term of the form yp y in Eq. ( 8.1 0| ) is opposite to the mixed 
term —ip v v in the differential operator (3.4) of the Fokker-Planck equation. This order is necessary to guarantee the 
conservation of probability. Indeed, multiplying the time evolution equation (8.9) by 5(y), and integrating both sides 
over x and y, the left-hand side vanishes. 

The correctness of th is or der can be verified by calculating the fluctuatio n det erminant of the path integral for the 
product of amplitudes ( |8.3| ) in the Lagrangian form, which looks just like (1.37), except that the difference between 
forward and backward trajectories y(t) = cc+(t) — %—(t) is nonzero at the endpoints. For the fluctuation which vanish 
at the end points, this is irrelevant. As explained before, the order is a short-time issue, and we can take tf, — t a — » oo. 
Moreover, since the order is independent of the potential, we ma y consider only the free case V( x ± y /2) = 0. The 
relevant fluctuation determinant was calculated in formula (4.5). In the Hamiltonian operator Q8.10p , this implies 
an additional energy —ij/2 with respect to the symmetrically ordered term r y{y,p y }/2, which brings it to ^yf> y , and 
thus the order in ( jg.lip . 



IX. CONCLUSION 



With the help of analytic regularization we have shown that the forward-backward path integral of a point particle 
in a thermal bath of harmonic oscillators yields, at large temperature, a probability distribution obeying a Fokkcr- 
Planck equation with the correct operator ordering which ensures probability conservation. By the same token, they 
yield the correct Langevin equations with and without inertia. 
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